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ABSTRACT 

The energetics of two numerical models of the ocean 
circulation are studied by use of available potential energy. 
Available potential energy, calculated by an approximation 
equation using a temperature perturbation field, compares 
favorably with that obtained using the exact expression 
using a height perturbation field. 

A general circulation model of the ocean is employed to 
study large-scale energy transformations, and a mid-ocean 
dynamics experiment (mesoscale) model is employed to study 
transformations of energy by stationary features and by 
transient motions. Results from the general circulation 
model are consistent with those of Haney (197^)» Results 
from the mesoscale model indicate that the transformation 
of energy is from available potential energy to kinetic 
energy of the transient eddy field. 
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I. INTRODUCTION 



A 'basic part of any study of an atmospheric or oceanic 
circulation system is a study of the behavior of the energy 
in the system and its transformation from one form of energy 
to another. The various forms of energy can be separated 
into three basic forms of energy; kinetic energy, potential 
energy and internal energy. Haurwitz (1941) showed that 
the potential energy and internal energy, in a vertical 
column, are proportional under hydrostatic equilibrium and 
thus may be combined into a single parameter, the total 
potential energy. Conversions of energy then occur between 
total potential energy and kinetic energy. 



A. AVAILABLE POTMTIAL ENERGY IK ATMOSPHERIC CIRCULATION 

Lorenz (1955)» in a study of atmospheric circulation, 
showed that since total potential energy was three to four 
orders of magnitude larger than kinetic energy, transformation 
of energy between the'tv/o energy forms v;as difficult to study. 
In a horizontally stratified atmosphere, heating of the 
atmosphere disturbs the stratification, allowing potential 
energy to be converted to kinetic energy. However, cooling 
of the atmosphere, which removes total potential energy from 
the system, also disturbs the stratification and creates 

conditions for conversion of total potential energy to kinetic 
energy. 
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Lorenz (1955) defined available potential energy as the 
difference between the total potential energy of the whole 
atmosphere and the total potential energy that would exist 
if the mass of the atmosphere were redistributed v/ith 
conservation of potential temperature to form a horizontally 
stable atmosphere. Then, under adiabatic flow, the sum of 
the available potential energy and kinetic energy is 
conserved. 

B. AVAILABLE POTENTIAL ENERGY IN OCEANIC CIRCULATION 

The use of the concept of available potential energy 
can be carried into the study of oceanic circulations. 

Several studies of the energetics of the ocean circulation 
systems have been based on available potential energy, as 
in Barrett (1971), and Wright (1972). The study by V/right 
dealt with changes in available potential energy in the 
ocean general circulation. The cooling of cold, high latitude 
water and the surface warming of warm equatorial water 
produces available potential energy by increasing the 
stratification of the ocean mass. Conversely, the sinking 
of the cold surface v;ater in the north tends to convert 
available potential energy to kinetic energy as it returns 
the water mass to a more stably stratified system. Barrett 
(1971) studied the available potential energy in the transient 
motion of Gulf Stream rings. Results of the study indicated 
that the decay of cold-core rings in the western Sargasso 
Sea region reduced the available potential energy of the 
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region. The presence of decaying cold-core rings in the 
western Sargasso Sea should have resulted in a rising of 
the mean thermocline in that region. Since the thermocline 
was not observed to rise significantly in the v;estern 
Sargasso Sea region, a local energy source must have been 
present to counteract the energy loss. 

C. AVAILABLE POTEI^TIAL ENERGY IN NUMERICAL OCEAN 

CIRCULATION MODELS 

This study applied available potential energy calculations 
to two experiments with a numerical model of the ocean general 
circulation; one simulating large-scale, quasi-steady state 
motion and the other simulating mesoscale transient motion. 

It v;as desired to calculate available potential energy 
using a perturbation height field compared with an approximate 
calculation using a temperature perturbation field. An 
expression for the balance of available potential energy was 
obtained from this approximation equation. The importance 
of the terms in the expression for the balance of available 
potential energy, and their corresponding physical processes, 
could then be examined in the two numerical experiments. 
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II. DESCRIFTION OF THE NUMERICAL MODELS 



The first experiment utilized the general circulation 
model of the ocean, as described in part A, for large-scale 
quasi-steady state motion. The Mid-Ocean Dynamics Experiment 
model (Robinson, et.al., 1974, MODS, Hot Line Nev/s , No. 57.), 
as described in part B, was utilized for the transient motion 
experiment. 

A. GENERAL CIRCULATION MODEL OF THE OCEAl^ 

The model represents an idealized ocean basin with a 
flat bottom and regula.r coast. It is everywhere 4 hin deep, 
90° longitude wide and extends from 5^.25S to 48.75N. 

The grid intervals are 3° of longitude and 2.5*^ of latitude. ' 
There are six vertical levels, placed at 20, 100, 300, 700, 

1500 , 3000 meters below the undisturbed surface. The model 

* 

is based on the primitive equations vrhich are integrated in 
time from specified initial conditions. The equations are 
based on the hydrostatic and Boussinesq approximations. 

The ocean is thus assumed incompressible. Salinity is 
neglected and density is a linear function of temperature. 

The equations are 




( 1 ) 
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( 2 ) 



dv~ _ -1 

d-L ' bf 
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( 5 ) 
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1 
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b (o-COSf)l r 0 

-I 


( 4 ) 


dr - 
dt 






-f- <^Ct) 


( 5 ) 



. e=e«(i- ”<cr-Tj) , (6) 



where the symbols have the following meanings; 

t = time 
'X = longitude 
^ = latitude 
Z = height 

A = angular speed of earth's rotation 
f = Coriolis force (=2^1 sin<^ ) 
a = radius of earth 
g = acceleration of gravity 
u = eastward velocity component 
V = northward velocity component 
w = vertical velocity component 
T = temperature 

T^= constant reference temperature 
Q = density 
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' -4..* 



p^= density of water at reference temperature 



oC = coefficient of thermal expansion 

H = depth of ocean basin 

convective adjustment of T if lapse 
rate becomes unstable. 

The lateral boundary conditions on heat and momentum 
at all four v/alls of the basin are those of perfect insulation 
and zero slip. The ocean floor is thermally insulated, but 
free slip is allowed. Due to the flat bottom, the vertical 
velocity is zero there. The circulation is driven and 
controlled by the wind stress and heat flux at the ocean 
surface. The wind stress is larger in the southern hemisphere 
and meridional wind stress is neglected'. The net dovmward 
heat flux at the ocean surface is expressed as 



atmospheric equilibrium temperature approximately equal to 
the air temperature at 10 meters and is a coupling 
coefficient determined from atmospheric climatological 
data. In this model, T is the temperature at the top level, 
20 meters below the surface. The model is a coarse-grid, 
high-viscosity model designed to show large scale ocean 
circulation. Complete documentation of the model is given 
in Haney (1974). 




* - -- 

where T^ is the ocean surface temperature, T^ is an 
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B. MID-0C3/iN DYNAMICS EXPERIMENT MODEL 



The mid-ocean dynamics experiment, MODE, model is a 

modification of the previously described model. It is 

designed to resolve mesoscale eddies with a horizontal 

grid resolution of 40 km and 5 levels in the vertical. A 

2000 km square section of the general circulation model area 

is represented by the MODE model, located along the western 

boundary and centered at 35N. Circulation is driven and 

controlled by the v/ind stress and thermal forcing, the same 

as in the previous model, but the eddy viscosity and 

conductivity were reduced to achieve mesoscale eddies. 

Eastern and v/estern boundaries are thermally insulated, 

no-slip walls, as in the general circulation model. Northern 

and southern boundary conditions are taken as free slip, 

zero normal velocity and a vertical temperature distribution 

which is independent of longitude and time. The ocean 

floor is thermally insulated as in the general circulation 

model. Initial conditions for the model were zero velocity 

and a temperature of 5^C at all levels. Twenty years of 

integration was performed on a 100 km X 100 km grid using 

7 2 

a horizontal viscosity of A^ = 2x10^ cm /sec and heat 

7 2 

diffusivity of Ajj = 1x10 cm /sec. Ten years of fine grid 

6 2 

calculations were then made using A = 2x10 cm /sec and 

m 

Ajj = 5x10^ cm^/sec. 
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III. DERIVATION OF AVAILABLE POTENTIAL MERGY EQUATIONS 



A. DOMAIN AVERAGED AVAILABLE POTENTIAL ENERGY 



The vertically averaged available potential energy 
contained within an idealized ocean basin may be calculated 
as follows. Assume the basin to be flat bottomed at a depth 
H. The total potential energy per unit volume in a column 
in the basin is given as 



P = 







where Z increases upward from the bottom where Z = -H. All 
variations in H are neglected. Neglecting also the dependence 
of density on pressure and salinity, the density may be 
written as a linear function of temperature, equation (6), 

The total potential energy per unit volume in the column 



is then given as 


£*o 




p r COtv)STAWT - 


1 f 


./ Z-drl. 


^ H 1 


L l-H V a / 





By defining Z(T) = 0 if T > T(0) and Z(T) = -H if T 1 T(-H) , 
the vertical average total potential energy per unit volume 
in the column may be vn?itten as _ 

- > •• T„i. a 

where Tmin and Tmax represent the coldest and warmest 

temperatures found in the basin respectively. The vertical 
average total potential energy averaged over the horizontal 
domain, denoted by ("") , is given by 
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P = COMSTAhJT 



H ^ 




2 

where the bar over Z denotes an average over an isothermal 

”2 

surface. The average potential energy is a minimum when Z 
is a minimum, ie. Z is equal to Z everywhere. Writing the 
height of a temperature surface as a mean height, Z, and 
a departure from the mean, Z*, then the domain average total 
potential energy is a minimum when Z' is zero everywhere. In 

“2 — ^ 

this case, Z can be replaced by Z . The minimum domain 
averaged total potential energy can then be written as 



- /■ 

Twim 

If the difference between the average bottom temperature, 
T(-H) , and the average bottom temperature in the state of 
minimum domain averaged total potential energy, T(-H)min, is 
assumed small and thus can be neglected, then the domain 
average available potential energy can be v/ritten as 




which can be written as 



A - 




( 7 ) 



Thus, the available potential energy is proportional to the 
variance of height of an isothermal surface. 
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B. APPROXIMTION TO THE DOMAIN AVERAGE AVAILABLE POTENTIAL 
ENERGY 

2 

The variance of height over an isothermal surface, Z’ , 
can be related to the variance of temperature over a horizontal 
surface, assuming the ocean is stably stratified. When the 
perturbation, Z*, from the mean height of an isothermal 
surface is greater than zero, the temperature of this region 
is less than the mean temperature by an amount T*. Similarly, 
when Z' is less than zero, T' is greater than zero. Therefore, 
changes due to Z' are partially compensated by changes in T*. 
The height of a temperature surface may be approximated by 

Z(T) Z(t) =£(t-t’) , 

which shows the height perturbation can be approximated by 

Z'(T) ^ . 

hr 

V/hen this is substituted into the expression for the domain 
average available potential energy, equation (7), the 
approximation to the domain average available potential 
energy is then 

which can be written as 

A 5 9<(Pp f */aT**_ dt . (8) 

H J S)T 

-M 

S2. 
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C. THE BALANCE OF AVAILABLE POTENTIAL ElIBRGY 



Since available potential energy is a very important 
source for kinetic energy, a study of the energetics of an 
ocean circulation should include a study of production and 
dissipation of available potential energy. Under adiabatic 
conditions an increase in kinetic energy is compensated by 
a decrease of available potential energy. When both energy 
forms increase, then formation of available potential energy 
must be occurring through diabatic processes. A study of 
these generation and dissipation processes and their 
importance may be made through use of the equation of the 
balance of available potential energy. The equation is 
derived as follows, based on the approximate form (8) of 
available potential energy. 

From equation (8), the time rate of change of available 



potential energy is 

,0 



- 8 «(e« f /li_ _ f . 

J)t ■ H HVSl/ J 






(9) 



By use of the thermodynamic equation (5)> equation (9) was 
rewritten to show the various terms responsible for the 
change of available potential energy. The term can 



be written as 



The term 



!(??) 



^(l. T'M ; t‘ Al’ . 

I at 

can be written as 



( 10 ) 









01 ) 
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The local rate of change of temperature and its average were 
obtained from the thermodynamic equation 

^ r _V (VT) t A.V'T 4 l( ^ + «T). oa) 

at ^2. “ 

The horizontal average of (12) gave 




where the terms - '^•(VT) and were zero since advection 
and diffusion of heat vanish at the boundaries of the domain. 
Equation (11) v/as then written using (13) as 



5 fc)T\ - _^r_ 5 (wt) 
atUz; ■ 2)21 Dz 




v/hich can be rev/ritten as 

A ( ^ f- 4 . 'r .1 

c>n iizi izl b'i e.<^J ’ 



where 






(14) 



(15) 



\ 



Equation (12) v;as v/ritten in terms of an average 
temperature and the departure from the average. Equation (13) 
was then subtracted from the result giving 

g’. -V(VT)-VCVf) - g)4(lt&)X^„V'T)' . 4«' ( . 6) 

Equation (.16) was then substituted into the right hand side 
of equation (10) to give 
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since 

T ^(wtV - o. 



The first term on the right side of (17) can be v/ritten as 



tT-v Cvt)- -tV vt - t'tvv -t'tS:*!, 

L SE J 81 

The term -TV*V T is equal to zero because V T is zero. By 
use of the continuity equation (if) , the equation reduces to 



tT- V - _ tV’ ^ . 

L 6-Z j ‘ 

Also, the second term on the right side of equation (1?) 
can be written as 



r-y*(\yT^) - 


T - . 


09) 


1 J 


^z\ i I 





The last term of equation (1?) can be written using ,the 
definition of Q as 




Substitution of equations (18), (19) > (20) into equation (1?) 
and using the resulting equation along with equation ( 15 ) 
results in the following form for equation (9) > the balance 
of available potential energy 
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(21) 




The balajice of available potential energy is due to 
five terms on the right side of equation (21). The first 
term is the transformation between available potential energy 
and kinetic energy due to the correlation between temperature 
and vertical velocity. The second term is a triple 
correlation term due to the convergence of the vertical 
transport of temperature perturbations. The third term is 
the generation of available potential energy due to a spatisil 
correlation between heating and perturbation temperature. 

The last two terms represent the change of available potential 
energy due to a change in the mean static stability. Lorenz 
( 1955 ) assumed the static stability v/as constant in time and 
also neglected the triple correlation term leaving only the 
first and third terms on the right side of equation (21). 

In the quasi geostrophic theory, the first and third terms 
are the only terms to appear in the balance of available 
potential energy. For the general circulation model, all 
five terms v/ere studied and for the mid-ocean dynamics 
experiment model the first three terms were studied. 
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IV. FINITE DIFFERENCE EQUATIONS 



The numerical program developed for calculations of 
available potential energy and the balance of available 
potential energy v/as designed for the general circulation 
model as developed by Haney (197^)* Therefore, the space 
finite differencing scheme and the boundary conditions 
were the sanie as those of the general circulation model, 

Haney (197^)» A staggered grid system was used for placement 
of the variables. Figure 1 shows the variable placement 
in the horizontal and figure 2 shows the variable placement 
in the vertical. Longitude positions are designated by 
the subscript i, latitude positions by the subscript j and 
vertical positions by k, increasing dov/nward. 

Horizontal averaging, denoted by q for some variable 
q. . v/as done in the following manner. Since the lateral 

J- J 

boundaries were located at integer i,j points, the mean 
value of a term to be averaged v;as calculated at the i+-^,j'‘'i- 
positions based on the terms value at, and the areas of, the 
grid points at the four surrounding i,j positions. These 

I 

mean values were then summed with equal weights to determine 
the domain horizontal average. 






LilTITUDE 




Figure 1 . Horizontal Placement of Variables 
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A. OTMERICAL CALCULATION OF AVAILABLE POTENTIAL ENERGY BY 
HEIGHT PERTURBATION 



To calculate the available potential energy by equation (7) 
it v/as necessary to calculate the temperature at the boundaries 
of the layers represented by the subscript k, represents 

the vertical mean temperature for a layer of thickness DZ^ 
where 




The temperature and height of the water column at the 

boundaries of each k level v/ere calculated such that T.. 

was maintained as a vertical mean* The actual temperatures 

^idk-|-* calculated using the following equations 

. 1 V r A ’ 3 I * 0 j i 

t-L — 

DZ^j, 






DZw ~ aKK • DZ KM-ya 
DH ’ 

where K14 represents the total number of levels in the model. 
The horizontal average of the vertical mean temperature, Tj, , 
was calculated as explained previously. The depth of T^^ 
was determined using 0.^,, ^ and Z, x use of the equation 

X J K— 2 



where is depth at v;hich the temperature 0 is 

equal to the horizontally averaged vertical mean temperature 
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in the layer k. The horizontal average of the depth of 
each of the mean temperature surfaces is cuid 

represents the average height of the isotherms betv/een 0^^ 
and 0j^^ ^ . The departure from the horizontal average was 
determined as the difference between the temperature surface 
and its horizontal average The resulting 

finite difference analog of (7) for available potential 
energy was of the form 




B. NUIvffiRICAL CALCULATION OF AVAILABLE POTENTIAL EI'IERGY BY 

TEl-lPERATURE PERTURBATION 

Equation (8) shov/s the analytical forms of the equation 
for the approximation of average available potential energy 
using the temperature perturbation, T*, instead of the height 
perturbation, Z', 

The horizontal average temperature, T, was calculated at 
each k level. By subtracting this value from the temperature 
at each grid position, a value for the temperature perturbation, 
T', was obtained at each grid position. Using these 

p 

perturbation temperatures, the term i-T* v/as calculated at 
each grid position and averaged to the half-integer points 
(i+i-, 0+i, k+i-). 

The term for the vertical derivative of the horizontal 
average temperature v;as calculated at the i,j integer points 
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as 



. . y r\J 



AX-k+y C ~ Tktl 

Z,.-Zv.. 

and averaged to the half-integer points. The approximation 



to the mean available potential energy was then obtained by 
combining these terms and summing over all levels as indicated 
by 






APE r y • 

H Z— \ 



(22) 



C. NUl-IERICAL CALCULATION OF THE BALANCE OF AVAILABLE 
POTEl^TIAL ENERGY 



The analytical form of the rate of change of available 
potential energy is shown by equation (9) as 



^ - 3^ 

■ 



/O- A-(-1,T \dY « ( . 

H ) J J 

“** iSz) •« [htl 



(9) 



The factors and in (9)j were calculated as in 

the numerator and denominator of (22). The remaining parts 
of equation (9) v/ere calculated by use of the right side 
of the thermodynamic equation as indicated by equations (10) 
and (11). Equation (10) was written to determine 

t>T’ 

by calculating T' — i at integer values of i,J,k and then 

a 

averaging to obtain its value at half-integer points. 
Equation (11) v/as calculated from 







c)t. ^ ^ 
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and then averaged to the half-integer grid points. The 
resulting finite difference form for equation (9)> the 
balcince of domain average available potential energy, is 



then 



. 3o(ec 



H I 




D. NUl-ffiRICAL CALCULATIONS OF THE TERiS FOR THE BALANCE 
OF AVAILABLE POTENTIAL ENERGY 

As indicated by equation (21), the balance of mean 
available potential energy contains five terms which are 
obtained by use of the thermodynamic equation (5)» 
Development of equation (21) showed that the five terms 
could be written as below 



-T W ^ = 


T* JJ- 


(18) 




T' ^-'7- (,\vt‘) - 


(19) 


C.c 




(20) 
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(24) 


ve,c I 


J* 


(25) 
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The terras in the brackets on the right side of the above 
five equations v;ere calculated using the identical finite 
difference method as used in the general circulation model 
to calculate the terras of the thermodynamic equation. The 
details are in Haney ( 1 974) • 
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V. RESULTS 



To nalie an accurate comparison of the two expressions 
for domain average available potential energy, equations (7) 
and (8), and to determine the importance of the terms 
occurring in the balance of available potential energy, 
equation (21), the general circulation model, as described 
in chapter II, was run for large-scale motion. A description 
of the conditions and the results obtained are given below. 

A. GENERAL CIRCULATION MODEL OF THE OCEAN 

To establish quasi-steady state conditions, the model, 

Haney (1974), was run in two phases. For the first phase, 
the linearized form of the vertical shear current was used. 

r 

The vertical shear current equations were obtained from the 
primitive equation form of the horizontal equations of motion, 
(1) and (2), by subtracting the vertical mean of equation (1) 
from equation (1) and by subtracting the vertical mean of 
equation (2) from equation (2). The vertical mean current is 
assumed constant in time. The temperature is determined from 
the complete thermal equation '5)« The first phase integration 
was for approximately 100 years. For the second phase, the 
complete nonlinear equations were used. The second phase 
integration v/as for 25 years. A more complete description of 
phase 1 and phase 2 may be found in Haney (1974). Twenty- four 
sets of model history data was taken at this point at intervals 
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of 25 days, representing a total of 600 days. The results 
showed the motion to be not strictly steady-state. As shown 
by Haney (197A-)» slight warming was occurring near the bottom 
after the integration period indicating a continuing downward 
flux of heat. Figure 3 shows the temperature at 20 meters 
below the surface, averaged over the last year of integration 
for phase 2. Figure 4 shows the vertical velocity at 60 meters 
below the surface averaged over the last year of the second 
phase. Comparing figures 3 s.nd 4 shows two large gyres 
centered on each side of the equator. In these gyres, warm 
water is being forced down. Near the equator and along the 
northern and southern boundaries, cold v/ater is being forced 
upward. Thus, available potential energy is being generated 




Figure 3. Temperature (°C) at 20 m, averaged over 
tne*last year of the second phase. 
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Figure 4, Vertical velocity (cm day" ) at 60 m, 
averaged over the last year of the 
second phase. Upward motion is 
shaded, 

1 , Results of Available Potential Energy Calculations 

Comparison of calculations of available potential 
energy as determined by equation (8) and the approximate 
form (9) showed that the available potential energy as 
calculated by the height perturbation to be slightly larger, 
(1383 ergs/cm^) , than the available potential energy as 
calculated by the temperature perturbation, (1232 ergs/cm^), 
figure 5* This discrepancy is perhaps due to Tmax and Tmin 
being not uniquely defined in Z coordinates. Fluctuations 
in available potential energy -due to a minor change in the 
model near the end of the 600 day period of record v/ere 
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Figure 5» APE by and by for the 

General Circulation Model. 
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simulated by both forms of available potential .energy 
calculations. Figure 5 indicates equation (8) calculated 
available potential energy to a sufficient degree of accuracy 
to be considered a good approximation equation for available 
potential energy in this model. 



2. Results of Calculations of the Balance of 
Available Potential Energy 



From equation (21), the balajice of available potential 
energy is represented by the sum of two integrals, 

= part 1 - part .2 



where 






1 . 





/ r yaT'~^ 1 rii_(wT ) ^ 

» J^l e,cJ 

4 . 5 . 






( 26 ) 



According to their definitions, Q' and Q contain the convective 
adjustment of temperature, ^ (T)', and S"«(l) respectively. 

In the calculations below, these convective adjustment terms 
were omitted due to the difficulty of recovering this 
irreversible process from the history records. 

Part 1 represents changes in available potential 
energy resulting from changes in the variance of the temperature 
perturbation on a horizontal surface. Heating of warmer 
regions and cooling of colder regions at the same depth 
increases the horizontal temperature variance and increases 
available potential energy. Part 2 represents changes in 
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available potential energy due to changes in mean static 

stability. Heating of upper levels and cooling of lower 

levels increases the static stability and decreases available 

potential energy. The magnitude of these integrals are 

shov/n by figure 6, Variations in these two parts after 

day 400 of the history record v/ere caused by the fact that 

the convective adjustment, ^ (T) , was changed to double 

c 

precision accuracy in the model. 

The processes represented by part 1 are serving to 

decrease the total available potential energy in this model, 

indicating the model, apart from ^ (T) ’ , is working to 

c 

reduce temperature perturbations. The magnitude of the 
three terms comprising part 1, equation (26), are shown 
in figure 7, where all three terms are of the same order of 
magnitude. The first term of part 1 shows the transformation 
of energy betvireen available potential energy and kinetic 
energy due to a correlation between the temperature 
perturbation and vertical velocity. Results indicate that 
in the domain average, energy is being transformed from 
kinetic energy to available potential energy. This is 
consistent with the results of Haney ( 1974 )> as discussed 
above, where a transformation of Icinetic energy of vertical 
shear motion to potential energy v/as occurring. Term 2 
represents a triple correlation due to convergence of the 

vertical flux of temperature fluctuations. Being a triple 
correlation of perturbation quantities, Lorenz ( 1955 ) 
assumed this term to be negligible in the atmosphere. In 
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of Eq* (26) vs. Time. 
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Figiire ?♦ Contribution by Terml(»), 
Term2(»), Term3(») to the 
Balance of APE, Eq. (26). 
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this general circulation model, term 2 is smaller than all 
of the other terms by only a factor of 2 to 3 'and thus it cannot 
be neglected. This term arises because the static stability, 
^T/SZj is not constant in Z, If it were, the term would 
vanish. The term would also vanish if the triple correlation 
were small. One reason it v/as not small is because of the 
forced nature of the vertical velocity, w, and temperature 
fields near the boundaries. With heat flux through the 
boundaries, neither the vertical velocity nor the temperature 
gradient in the horizontal v/ould have been as large as 
occurred. The third term of part 1 shows generation of 
available potential energy due to a spatial correlation 
between heating and perturbation temperature. This term 
T/as calculated without the contribution from $ (T) ' , 
present in equation (20). The negative result of term 3 
indicates a loss of available potential energy due to the 
non-convective heating. The magnitude of this term indicates 
the importance of the thermal diffusion terms in this 
oceanic general circulation model. 

The processes represented by the terms of part 2, 
figure 8, were of the same order of magnitude, and tend to 
increase the available potential energy in the model. The 
first term of part 2, term 4, is due to the vertical flux 
of mean temperature and the second term, term 5> is due to 

mean heating. As in term 3 » the convective adjustment 

term, S (t)', v;as absent from the mean heating term. ^ '(t) 
c c 

tends to increase the mean static stability, ^T/^Z, and to 
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Balance of APE, Eq. (26). 
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decrease available potenticil energy. It's inclusion would 
allow terms Zf and 5 to sum to a negligible value. Lorenz 
(1955) assumed part 2 was zero for steady-state atmospheric 
motion. Under conditions of seasonal radiation heating, 
the vertical flux of mean temperature, term may not be 
in phase with the mean heating that occurs, term ajid 
thus v/ould not cancel. The general circulation- model was 
not designed to sho\7 seasonal radiation heating. 

The energy flow diagram depicting terms 1 - 5 in 
equation (26) is shovm in figure 9» Terms l,Zf and 5 were 
serving to increase available potential energy by an amount 
of 1.705x10"^ ergs/cm^sec. Terms 2 and 3 were serving to 
decrease available potential energy by an amount of .908x10”^ 
ergs/cm^sec. The available potential energy during the 
record period v/as observed to be decreasing, figure 5j with 

■z 

an average value of 1232. ergs/cm*^. The discrepancy between 
the terms serving to increase available potential energy, 
part 1 - part 2>0, and the observed decrease in available 
potential energy, figure 3> was determined to be due to not 
including the ^ (T) term in the balance of available potential 
energy calculations. The value of the term required to 
bal£ince the energy flow diagram, as in steady-state conditions, 
is indicated by term 6 and shovm by the dashed line in 
figure 9 to be approximately .797x10"^ ergs/cm^sec. These 
results indicated convective adjustment was an important 
process in the model, as was expected since the model 
contained a large region of near isothermal water at high 




Figure 9* Energy Flov/ Diagram for the 
General Circulation Model. 



43 



latitudes where convective adjustment would frequently 
occur. .The convectj.ve overturning plays no role in the 
kinetic energy budget and therefore the arrow for term 6 
shov;s this effect to be an external loss to the system. 

B. MID-OCEAN DYITAI-IICS EXPERIMENT MODEL 

Figure 10 below shows the instantaneous streamfunction 
for day 10080 during the seventh year of fine grid 
calculations. Tv;enty-four sets of data for the experiment 
were taken from day IOO 5 O to 10110 at 2x day intervals. 




Figure 10. Instantaneous streamfunction for 
day 10080 of the Mid-Ocean 
Dynamics Experiment Model. 
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To determine the available potential energy as a result 

of transient motion and standing motion, equation (8) was 

written with temperature in terms of a time average 

temperature, CT], for standing motion, and a departure from 

* 

the time average, T , for transient motion. Since available 
potential energy is determined from the temperature 
perturbation, the horizontal average of the time averaged 
temperature v/as determined. Then the perturbation temperature 

* I 

of the standing motion, \^Tj * and transient motion, T , was 
determined as below 

T = t t” 

T = T' •* tT]' f T* J (27) 

where lT3 is the horizontal averaged standing temperature 
and T is the horizontal averaged transient temperature. 

The term J ’ is the standing eddy temperature and T is 
the transient eddy temperature. The eddy temperature is 
then the sum of the last two terms of equation (27). 
Substitution of the eddy temperatures into equation (8) 
resulted in the following form for domain average available 
potential energy. 
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The time average4 available potential energy was written 
as 



H t a. [-] 



) 



( 29 ) 



p 'K' i-| 

where the term \|-1 lT" 1'T ] vanishes. In this model, the 

static stability factor T/^Z} was assumed constant in 
time and written as If . Equation (29) shows time average 
available potential energy as the sum of a part due to the 
standing field and a part due to the transient motion. 

For the study of the balance of available potential 
energy in the mid-ocean dynamics experiment model, only 
the first three terms of equation (21) were considered 
important for further study. Because lateral boundary 
conditions allowed heat to flow in or out of the model, it 
was reasoned that convective adjustment would not play an 
important role in the model. With little convective 
adjustment and steady forcing, it was therefore assumed 
that the long-time average of the sum of terms if and 
part 2 of equation (26), would be negligible. Calculations 
of the balance of total available potential energy, belov/, 
showed the terms of part 2 summed to approximately zero 
for each set of data as v;ell as in a time average so the 
assumption was Justified. The study therefore concentrated 
on transformations of energy due to -JT* . From the first 
three terms of equation (21), the balance of available 
potential energy was approximated as 
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( 50 ) 






f-Tw - > ^ tv*) 

..'.t (g)»U J (||)e.cj 

1 . 2 . 5 . 






Equation (30) was written in terms of a part due to transient 
motion, q , and a part due to standing motion, M. for a 



quantity q and then time averaged as 



raAl , 3He./'f-[wT[Tr 

liTj ' => 



la. 



lb. 



2a. 



2b. 



Si / y ^ 

2c. 2d. 33-» 5b. 

where the quantities containing only one term due to transient 
motion were zero in the long-time average. The values of 
terms 2c and 2d were not calculated in this study but their 
sum could be inferred from terms 2, 2a and 2b. Since heat 
diffusivity v/as allowed across the boundaries, the definition 
in (25) was changed to the following 

1. ICb 



h ■ • (" S) 



since AjjV^T was not equal to zero in the MODE model. 



1 . Results of Available Potential Energy Calculations 



Figure 1 1 shov/s the total available potential energy 
as determined from equation (28) to be slowly increasing 
during the history record. The first term on the right side 





Figure 11. APE vs. Time for the MODE Model 
(Total Motion) 
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of equation ( 28 ), due to standing motion was constant in time 
and had a value equal to 116»65 ergs/cm^, shown in figure 11. 
This v/as an order of magnitude smaller than the available 
potential energy of the general circulation model and it 
points out the importance of the horizontal diffusion of 
heat across the northern and southern boundaries in preventing 
large horizontal gradients of temperature. The. third term 
of equation ( 28 ), due to transient motion, is shovm in 
figure 12. This term appears to have a period of approximately 
100 days, but is too small in amplitude to significantly 
affect the total available potential energy. Since the left 
side of equation (28) v/as known from figure 11, and two 
terms on the right side of equation ( 28 ) v/ere known, the 
middle term on the right side of equation (28) was found to 
be primarily responsible for the slow increase of available 
potential energy seen in figure 11. The time average of 
the total available potential energy, corresponding to the 
left side of equation (29) was 117.3 ergs/cm^. The first 
term on the right side of equation (29), the available 
potential energy due to stauiding motion, is the same as the 
first term on the right side of equation ( 28 ), equal, to 
116.65 ergs/cm^. Thus the time average available potential 
energy due to trauisient motion, the last term of equation (29), 
had a value of .66 ergs/cra^ and was only a small part of the 
time average total available potential energy. Being so 
small, the available potential energy due to transient 
motion would probably not be noticed on a temperature map 
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Figure 12. APE by iT‘^ vs. Time for the MODS Model 

(Transient Motion) 
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from the model. If the variations in kinetic energy were of 
this order of magnitude, it could be important in the ld.netic 
energy balance which v;ould be subject to future investigation, 

2, Results of the Balance of Available Potential 
Energy Calculations 

By the use of flow diagrams, the magnitude of the 
terms responsible for the balance of available potential 
energy can be easily seen. Figure 13 shows the balance of 
the time average total available potential energy, 
corresponding to equation (30), This figure indicates energy 
was being transformed from kinetic energy to available 
potential energy by term 1 of equation (30) at a rate faster 
than it was being removed by terms 2 and 3» The average 

rate of increase of total available potential energy, as 

•“5 3 

shown in the diagram, was approximately ,0if6xl0 ergs/cm'"^sec. 

From the observed change in available potential energy, 

figure 11, there was an average rate of available potential 

energy change during the 60 days of ,05x10 ergs/ciar sec. 

Figure Vi shows the value of the energy transformations 

due to standing motion, corresponding to terms la, 2a, 3a of 

equation (31 )> and figure 15 shov/s the value of the energy 

transformations due to transient motion, corresponding to 

terms 1b, 2b, J>b of equation (31 )• The average rate of 

increase of available potential energy by the standing field, 

as shown in the diagram, was approximately ,075x10 

ergs/cm^ sec. The average rate of decrease by the transient 

-5 / 3 

field, as shown in the diagram, was ,035x10 ergs/cvr sec. 
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Figure 13» Energy Flow Diagram for the MODS Model 
(Total Motion Eq. 30) 






sec 

sec 
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Net Change of APE +0. 075x1 ©"’^ergs/cm^sec 



Figure 14. Energy Flow Diagram for the MODE Model 
(Standing Motion, Eq. 31) 
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Figure Energy Flow Diagram, for the MODE Model 
(Transient Motion j Eq. 50 
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Instead of such changes in the available potential energy 
fields, this perhaps simply indicates that a conversion of 
available potential energy from the standing to the transient 
fields was occurring. The conversion term for available 
potential energy between the standing and transient fields 
was not calculated because it also contained terms of 
standing and transient motion as doss terms 2c and 2d. The 
conversion terra must have had a value between .055 to .075 
xlO ^ ergs/cm-^ sec. Summing of the net increases by the 
standing field and the transient field gave a net increase 
for the total of .0^0x10**^ ergs/cm^ sec, v;hich compares 
with the value of .046x10*"^ ergs/cm^ sec seen in figure 15 . 

The sura of terras 5a, figure IZfj and 5t>, figure 15> was equal 
to term 5 of figure 15 as it should be, since in the time 
average term 5 is composed of only term 5a and 5t. Similarly, 
term 2a, figure 14 > and 2b, figure 15 > sum to equal term 2 
of figure 15. However, from equation (51 ) term 2 is composed 
of four terms in the time average. Therefore, terms 2c and 
2d must be negligible or sum to a negligible value. 

The sum of la and lb is less than term 1 by an 
amount of .006x10 ^ evgs/ cnr sec. Some round-off errors 
may have occurred to cause the small discrepancy. Comparing 
figure 15 and figure 14 , the total energy being transformed 
from IdLnetic energy to available potential energy, term 1 , 

was almost totally by the standing field, terra la. Since 
terra lb, figure 15 , v;as small and positive, the source of 
kinetic energy for the transient motion was not available 
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potential energy during the 60 day period of the MODS model. 
Figure- 16 shows the total transformation from kinetic energy 
to available potential energy, term 1 , v/as undergoing large 
fluctuations. At some times the transient eddies did draw 
energy from available potential energy, since 1b, figure 17, 
v/as at times negative, but the energy flow in either direction 
was small. These results may indicate that the- transient 
eddies were not baroclinic for the 60 day period. The 
amplitude of the eddies of the MODS model v/ere known to be 
large. With large amplitude eddies, the dissipative forces 
were expected to be large, indicating the transformation of 
energy from the source of energy for the Icinetic energy of 
the eddies v/as large. Therefore the primary source of energy 
for the transient eddies was not available potential energy 
during the 60 day history record. 

It was suspected that errors may have occurred due 
to the S-g- day rate of sampling the data. Figure 16 shows the 
three terms of equation (30) as a function of time. Term 1 
and 2 fluctuate to a greater extent than term 3» Since 
both term 1 and 2 are correlated v/ith vertical velocity, w, 
the possibility of aliasing is further supported. Fluctuations 
in term lb, figure 17, and 2b, figure 18, also lends support 
to this possibility due to their fluctuations. Neither 
term 3 of figure 16 nor term 3D of figure 19 showed a large 
degree of fluctuation and neither of these terms were 
correlated with vertical velocity. Therefore, the above 
energetics v/ere recalculated using a 120 day history record 
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Figure 16. Contribution by TermlxlO^(») , 
Term£xlO^(«) , Terra3xlO^(#) , to the 
Balance of APE, Eq, (30) • 
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Figure 1?. Contribution by Termlb xlO^ 
to the Balance of APE. 



58 



Vosinij 




Figure 18. Contribution by Term2b xio"^ 
to the Balance of APE. 
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v/hich began on day 100^0 and ended on day 10170. Ko essential 
changes occurred in any of the terms except term lb, which 
was still small, but negative for the 120 day average. Since 
a 120 day average would probable be more accurate than a 60 
day average, the eddies were felt to be baroclinic. 
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VI. CONCLUSION 



The results obtained from this study indicated that 
both the general circulation model and the MODE model were 
mechanically driven by the wind stress and not thermally 
driven to a great extent. The below diagram shov/s the 
transformation of energy from its generation by the v/ind 
stress to its dissipation through transient eddies as 
indicated by the results from the MODE model. 




The atmospheric wind circulation acts on the ocean surface 
creating large scale standing motion. Some of the kinetic 
energy of the standing motion is transformed to standing 
available, potential energy of the ocean. Approximately 
half of the input of energy to standing available potential 
energy is lost through non-convcctive heating by standing 
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motion. About one-fourth of standing available potential 
energy is transformed to transient available potential energy, 
where a large portion is lost through transient non-convective 
heating. A small portion, in the transient available 
potential energy balance, is transformed to kinetic energy 
of transient motion v/here it is lost through dissipation. 
Whether the transient available potential energy is the 
primary source for the transient kinetic energy must be 
determined by a study of the transient kinetic energy 
balance. 

In the atmosphere, forcing is by the heating and the 
energy transformations are from available potential energy 
of standing motion to transient available potential energy, 
then to transient kinetic energy. 

The main conclusion is that baroclinic instability is 
a possibility for maintaining the MODE eddies. The general 
circulation model showed energy transformations from eddy 
kinetic energy to eddy available potential energy for large- 
scale quasi-steady state motion. The MODS model shov/ed the 
energy transformations to be opposite, from eddy available 
potential energy to eddy kinetic energy, for transient 
motion. Therefore, during periods of baroclinic instability, 
eddy available potential energy is aiding in the maintenance 
of the eddies. 
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